Direct determination of the strange and light quark condensates from full lattice QCD 
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We determine the strange quark condensate from lattice QCD for the first time and compare its 
value to that of the light quark and chiral condensates. The results come from a direct calculation of 
the expectation value of the trace of the quark propagator followed by subtraction of the appropriate 
perturbative contribution, derived here, to convert the non-normal-ordered mipip to the MS scheme 
at a fixed scale. This is then a well-defined physical 'nonperturbative' condensate that can be used 
in the Operator Product Expansion of current-current correlators. The perturbative subtraction is 
calculated through 0(a s ) and estimates of higher order terms are included through fitting results at 
multiple lattice spacing values. The gluon field configurations used are 'second generation' ensembles 
from the MILC collaboration that include 2 + 1 + 1 flavors of sea quarks implemented with the Highly 
Improved Staggered Quark action and including u/d sea quarks down to physical masses. Our results 
are : (ss) W (2GeV) = -(290(15) MeV) 3 , (Zi) W (2GeV) = -(283(2) MeV) 3 , where lis a light quark 
with mass equal to the average of the u and d quarks. The strange to light quark condensate ratio 
is 1.08(16). The light quark condensate is significantly larger than the chiral condensate in line with 
expectations from chiral analyses. We discuss the implications of these results for other calculations. 



I. INTRODUCTION 

A critical feature of the nonperturbative dynamics of 
QCD at zero temperature is the condensation of quark- 
antiquark pairs in the vacuum, spontaneously breaking 
the chiral symmetry of the action. The value of the chiral 
condensate (the quark condensate at zero quark mass) is 
then an important parameter for low energy QCD Q]. 
The well-known Gell-Mann, Oakes, Renner (GMOR) re- 
lation [2]: 

JlMl = m u + m d (0\uu + dd\0) 

4 2 2 U 

connects the u/d quark masses times condensate to 
the square of the mass times decay constant for the 
Goldstone boson of the spontaneously broken symmetry. 
Eq.([T]) has normalisation such that f n = 130 MeV. The 
GMOR relation holds in the limit of m u ,m,i — > 0. A 
value for this chiral condensate can be derived from the 
chiral extrapolation of lattice QCD results for light me- 
son masses and decay constants. See, for example, the 
recent result of —(272(2) MeV) 3 for the chiral condensate 
in the MS scheme at 2 GeV using SU(2) chiral pertur- 
bation theory [3]. 

The determination of the quark condensate for non- 
zero quark masses is more problematic because, depend- 
ing on the method used, there are various sources of un- 
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physical quark mass dependence and a careful definition 
of the condensate is required. This definition must be 
phrased in terms of the Operator Product Expansion 
(OPE) since this is the context in which the conden- 
sate appears [U HI |5] . The OPE allows us to separate 
short and long-distance contributions in, for example, a 
short-distance current-current correlator. The expansion 
is in terms of a set of matrix elements of local opera- 
tors multiplied by coefficient functions. The aim is for 
all the long-distance contributions (with scale < fi) to 
be contained in the matrix elements and the short dis- 
tance contributions (with scale > /i) in the coefficient 
functions. A key matrix element, since it corresponds to 
a relatively low-dimensional (d = 3) operator, is that of 
the quark condensate. The clean separation of scales in 
the OPE only works if the local operators are not normal 
ordered [BJ |7]- Then the coefficient functions are ana- 
lytic in the quark masses and therefore free of infrared 
sensitivity. This means, however, that the quark mass 
dependent mixing of rmpip with the unit operator must 
be taken into account and that the vacuum matrix ele- 
ment of rmpip is not cut-off independent. The quantity 
that appears in the OPE is the vacuum matrix element 
in, for example, the MS scheme at the scale /x. We can 
derive this matrix element from lattice QCD and we give 
results here for /i = 2 GeV. The results can easily be run 
to other scales, as appropriate. 

The value of the condensate for quarks of non-zero 
mass up to that of the strange quark is needed in a num- 
ber of calculations involving light quark correlators. In 
lattice QCD it is frequently easier and statistically more 
precise to use strange quarks than very light quarks in 
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contexts where the quark mass is not expected to be im- 
portant. Then the strange quark condensate appears in 
the calculation, however. Examples include the matching 
to continuum QCD perturbation theory of lattice QCD 
calculations of moments of heavy-light meson correla- 
tors [8] and of light meson correlators at large space-like 
q 2 [5] • Such calculations are used to extract quark masses 
and the strong coupling constant, a s . A continuum ex- 
ample where the strange quark condensate is needed is in 
the determination of the strange quark mass, m Sj from 
hadronic r decays [TU] , 

Current estimates of the value of the strange quark 
condensate vary by almost a factor of two [HJ [12] . It is 
not even clear whether the strange condensate is larger 
or smaller than the light quark condensate. For very 
large quark masses, m q > Aqcd, say, so that the quark 
mass dominates the propagator, it seems clear that the 
condensate should fall to zero, but this does not help 
in determining the slope of the condensate with m q for 
small quark masses. 

Here we address the determination of the strange (or 
other non-zero mass) quark condensate by direct calcu- 
lation in full lattice QCD. By direct we mean that we 
determine the vacuum expectation value of the strange 
quark propagator as well as the light quark propagator 
on a range of gluon field configurations at different values 
of the lattice spacing and sea quark masses. To isolate 
the low-energy nonperturbative value of the condensate 
from these results requires the subtraction of a perturba- 
tive contribution. The perturbative contribution in lat- 
tice QCD has two pieces. One diverges as a — > and 
dominates the vacuum expectation value of the strange 
quark propagator, particularly on our finer lattices. The 
second piece contains infrared sensitive logarithms of the 
quark mass which cancel against similar terms in contin- 
uum perturbation theory allowing an infrared safe defi- 
nition of the condensate for use in the OPE, as discussed 
above. 

The error in the final result then depends on how well 
this subtraction can be done. Here we use an explicit 
calculation of the perturbative pieces through 0(a s ) and 
fit for unknown higher order terms. The known quark 
mass and a dependence of these unknown terms helps 
in constraining them along with the very small statisti- 
cal errors in our lattice results. We also use a particu- 
larly good discretisation of the Dirac action known as the 
Highly Improved Staggered Quark (HISQ) formalism [T3] 
on 'second generation' gluon field configurations so that 
discretisation errors in the physical nonperturbative re- 
sults are small. 

The paper is laid out as follows. In Section [H] we de- 
scribe the theoretical background to direct calculations 
of the quark condensate in lattice QCD. Section [ITT] gives 
our lattice QCD results on gluon configurations with 
2+1+1 flavors of sea quarks, describing the calculation 
of the perturbative contribution that is subtracted and 
then the procedure for fitting the remaining nonperturba- 
tive condensates as a function of quark mass and lattice 



spacing. We also give results from configurations includ- 
ing 2+1 flavors of sea quarks over a wider range of lattice 
spacing values but studying only the strange quark con- 
densate in Appendix[Dj In Scction [TV] we compare to pre- 
vious values and discuss the implications of our results 
for both zero and finite temperature QCD calculations. 
Section [V] gives our conclusions. 

II. THEORETICAL BACKGROUND 

The direct determination of the chiral condensate in 
lattice QCD requires the calculation of the expectation 
value over an ensemble of gluon fields, U, of TrM -1 
where M is the lattice discretisation of the Dirac matrix. 
The quark action for a given quark flavor, 

S f = i>M f ^ (2) 

and 

$V»> = (0|^/|0) = -^(TvMfiuy^u, (3) 

where the trace is over spin, color and space-time point 
and the gluon fields in the ensemble used for the average 
include the effect of sea quarks (of all flavors, not just /) 
in their probability distribution. V is the lattice volume, 
L z x T. For a naive discretisation of the Dirac action M 
takes the form: 

M = 7 M A M + m (4) 
where is a covariant finite difference on the lattice: 

and m is the quark mass for that flavor. Because of 
fermion doubling this formalism describes 16 'tastes' of 
quarks in 4-dimensions rather than just 1 and we must 
divide the right-hand side of Eq. ^ by the number of 
tastes, N t = 16. The staggered formalism is derived from 
this naive formalism by a rotation which allows the spin 
degree of freedom to be dropped. In that case the quark 
field becomes a 1-componcnt spinor, which is numerically 
very efficient, and Nt = 4. 

For m = the eigenvalues of M for either naive or 
staggered quarks, are purely imaginary and come in ± 
pairs. Therefore, in the absence of exact zero modes, 




A calculation at m = on a finite volume lattice would 
then give an answer for the quark condensate of zero. 
This does not mean that chiral symmetry is unbroken, 



3 



however. The problem arises because the broad distribu- 
tion of non-zero eigenvalues (ignoring topological near- 
zero modes) drops to zero near the origin in a way that 
depends strongly on the volume. Once the quark mass is 
below the minimum of the non-zero eigenvalues the result 
for (ipip) wm be distorted. A more careful consideration 
of limits must be made. If V is taken to infinity before m 
is taken to zero then the sum over eigenvalues can be re- 
placed by an integral and the Banks-Casher relation [2] 
is obtained. This connects the zero-mass condensate to 
the spectral density at the origin: 



£ = -(ipip(m = 0)) = 



MO) 



(7) 



Thus S can be obtained from studies of the spectral 
density and, more recently, has also been obtained from 
matching the distribution of low eigenmodes to random 
matrix theory in the e regime |15H17j . Here we are 
more concerned with extracting a condensate at non- 
zero quark mass, for example at the strange quark mass, 
and so the issue above is not relevant. We will work 
on large volume lattices (over ten times larger than the 
study in [TS] that looked at staggered eigenvalues in 
the e regime) at quark mass values that are well within 
the distribution of non-zero eigenvalues. Our results for 
(TrM -1 ) then include both the effects of a non-zero value 
for p(0) and a non-zero quark mass. 

A second issue arises, however, in the extraction of 
a physical nonperturbative value for the condensate at 
non-zero quark mass. A perturbative contribution ap- 
pears from mixing between the scalar iptp operator and 
the identity since the identity operator has a vacuum ex- 
pectation value in the trivial perturbative vacuum. This 
perturbative contribution vanishes at zero quark mass 
since chiral symmetry is not broken in perturbation the- 
ory (for the same reason as that given on the lattice above 
in Eq. Q). At non-zero quark mass it contains odd pow- 
ers of to starting with a quadratically ultra-violet diver- 
gent term linear in to. This can be illustrated with a 
tree-level calculation in the continuum to give: 



(iH>) = 



d*k 12m 
o (2tt) 4 k 2 + m 2 



(8) 



— mA z + til log — — 
4tt 2 V A 2 + 



The quadratic ultraviolet divergence depends on the 
scheme used but the to 3 log(m/A) term is universal since 
it arises from the infrared part of the integral. The coeffi- 
cient above agrees with that obtained for the MS scheme 
in [THj and on the lattice for highly improved staggered 
quarks to be described in section |III A| We stress that 
a perturbative contribution of this kind is present for 
all lattice regularisations of QCD, whatever their chiral 
symmetry properties, and so must be calculated and sub- 
tracted to give a physical result. The quadratic diver- 
gence present for Ginsparg- Wilson fermions is demon- 
strated in the quenched approximation in, for exam- 



ple, [2D] and the additional divergences for the Wilson 
formalism with broken chiral symmetry in |21j . 

This subtraction is somewhat analogous to subtracting 
perturbative contributions to the mean plaquette to ob- 
tain the nonperturbative gluon condensate. That, how- 
ever, is extremely difficult to do because the nonper- 
turbative condensate contribution to the plaquette is so 
small. This contribution is given at leading order by: 



SPcond = {a s G 2 /ir). 
3d 



(9) 



If we take the value of the gluon condensate as 0(Aq CD ) 
then (a s G 2 /n) sa 0.005 GeV 4 . On very coarse lattices for 
which awl GeV -1 , this contributes less than 1% to the 
value of the plaquette. On finer lattices the nonpertur- 
bative condensate contribution is even smaller because it 
falls as a 4 while the perturbative contribution falls only 
as a s (d/a) for some scale d. This means that the pla- 
quette is in fact a very good variable to use for the de- 
termination of a s from lattice QCD calculations but not 
for the determination of the gluon condensate [22] . For 
larger Wilson loops the gluon condensate contribution is 
larger, being proportional to the square of the area of the 
loop, but the coefficients in the perturbative series also 
become larger. 

The determination of the nonperturbative quark con- 
densate from (TrA/ -1 ) is in much better shape than this 
for several reasons. The main one is that the nonper- 
turbative condensate contribution to (TrM^ 1 ) in lattice 
units is only a factor of a 2 smaller than the leading 
perturbative contribution rather than a 4 . In addition 
the perturbative contribution is suppressed by the quark 
mass, which is small for the u/d and s quarks we will 
consider here. The perturbative contribution is a well- 
defined function of the quark mass at every order in per- 
turbation theory and so results at several values of the 
quark mass, and the lattice spacing, can be used to con- 
strain unknown higher orders, beyond the 0{a s ) that we 
have explicitly calculated, and we will make use of that 
here. 



III. LATTICE QCD CALCULATION ON 
n/=2 + l + l GLUON CONFIGURATIONS 

The gluon field configurations used here are listed in 
Table [T] They were generated by the MILC collabora- 
tion |23j using a tadpole-improved Liischer-Weisz gauge 
action with coefficients corrected perturbatively through 
0{a s ) including pieces proportional to n/, the number 
of quark flavors in the sea 24J. The gauge action is 
then improved completely through <D(a s a 2 ). Sea quarks 
are included with the highly improved staggered quark 
(HISQ) action [13] which has been designed to have very 
small discretisation errors. Discretisation errors are for- 
mally removed through 0(a 2 ) but higher order errors, 
particularly staggered taste-changing errors, are seen to 
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TABLE I: Details of the MILC gluon field ensembles used in 
this paper. /3 = 10 /g 2 is the 5*17(3) gauge coupling and L/a 
and T/a are the number of lattice spacings in the space and 
time directions for each lattice. ami iSen , am a , aea and am c , SM 
are the light (up and down taken to have the same mass), 
strange and charm sea quark masses in lattice units, a is the 
lattice spacing in fm determined from the decay constant of 
the r/ s meson in [25] with values for 3, 6 and 8 added here. 
The ensembles 1, 2 and 3 will be referred to in the text as 
"very coarse", 4, 5 and 6 as "coarse" and 7 and 8 as "fine." 



Set 


P a/fm 


am Uaea 


CLTYl s ^sea 




L/a x T/a 


1 


5.80 0.1546(11) 


0.013 


0.065 


0.838 


16x48 


2 


5.80 0.1526(8) 


0.0064 


0.064 


0.828 


24x48 


3 


5.80 0.1511(8) 


0.00235 


0.0647 


0.831 


32x48 


4 


6.00 0.1234(8) 


0.0102 


0.0509 


0.635 


24x64 


5 


6.00 0.1218(6) 


0.00507 


0.0507 


0.628 


32x64 


6 


6.00 0.1206(6) 


0.00184 


0.0507 


0.628 


48x64 


7 


6.30 0.0899(7) 


0.0074 


0.0370 


0.440 


32x96 


8 


6.30 0.0875(7) 


0.0012 


0.0363 


0.432 


64x96 



be smaller with HISQ than with the earlier asqtad stag- 
gered quark action [T31 [23] • The HISQ action used here 
has two smearing steps for the gluon field appearing in 
the quark action with a U(3) projection of the smeared 
links between the two smearing steps. The configura- 
tions include a sea charm quark in addition to up, down 
and strange. These configurations are then said to have 
2+1+1 flavors in the sea, since the u and d quarks are 
taken to have the same mass (denoted mi here). This 
is heavier than the average u/d mass in the real world 
on most of the configuration sets but there are three for 
which the u/d mass has its physical value (3, 6 and 8). 
The s and c masses are tuned as closely as possible to 
their correct values on each set. The tuning of the sea s 
quark mass is accurately done - typically to better than 
5% - so the u/d quark mass can be accurately calibrated 
in terms of the s quark mass for chiral extrapolations. 
The values of the lattice spacing for most ensembles were 
determined in [25] using the decay constant of the rj s me- 
son. The values vary from 0.15 fm to 0.09 fm as we go 
from the very coarse to the fine lattices. The spatial 
volumes are large, from (2.5 fm) 3 when mi/m s « 0.2 to 
(3.7fm) 3 when mi/m s « 0.1. 

On each of these ensembles we determine (TrM -1 ) for 
HISQ valence quarks for various quark masses. To do 
this we use an identity that relates the quark propagator 
for staggered quarks to a product of quark propagators: 

^TrAf^ 1 = ^[M^M-^i-ir 

% n 

= ^TrlikCl 2 (10) 



Here and n are arbitrary lattice sites and am q is the 
quark mass in lattice units used for the quark propagator. 
The righthand side of Eq. ( 10 ) is simply the correlator 




FIG. 1: Autocorrelation function, Cat of Eq. |12| for the 
strange and light quark condensates on coarse set 6 with phys- 
ical mass light sea quarks. The a;-axis, AT, is the separation 
in time units between configurations. The strange conden- 
sate results are given as blue crosses and the light condensate 
with red pluses. Errors in Cat are estimated by dividing the 
configuration time series into five consecutive sets. 



made of a quark and antiquark of mass am q . Summing 
over n projects on to zero spatial momentum and sums 
over timeslices. Thus, dividing both sides by 4, the num- 
ber of tastes for staggered quarks, we obtain: 



i 



(11) 



The raw condensate value on the left-hand side of this 
expression is normalised to the single flavor case and the 
pion correlator on the right-hand side is the usual zero- 
momentum Goldstone meson correlator. This allows us 
to determine by summing over the Goldstone pseu- 

doscalar correlators calculated in [25]. Eq. ( 10 1 is derived 
in [26] for a single propagator origin, 0, but the deriva- 
tion can trivially be extended to hold for the random 
wall source that we use for our correlators in [25j . The 
identity holds configuration by configuration for lattice 
QCD quark formalisms with sufficient chiral symmetry 
and in the continuum for a specific gauge field back- 
ground. We give an explicit proof of this in Appendix [A] 
Since our Goldstone pseudoscalar correlators are sums of 
positive numbers they are particularly precise and this 
precision then carries over to our condensate results. For 
the light condensate we use the pion correlators made of 
light quarks and for the strange condensate we use the rj s 
correlators made of strange quarks. We stress that what 



we calculate using Eq. ( 10 ) is the vacuum expectation 
value of the condensate (and not a specific 'in-meson' 
value) despite the fact that we determine it for conve- 
nience from a meson correlator. 

Table |H] gives the valence quark masses used in our cal- 
culation and the raw results for the condensate obtained 
from Eq. (10 1. The correlator calculations used 16 'ran- 



between and n for the Goldstone pseudoscalar meson dom wall' time sources on approximately one thousand 
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TABLE II: Raw (unsubtracted) values for the light and strange quark condensates in lattice units calculated for valence masses 
given in columns 2 and 3. The results use the correlators calculated in [25] (via Eq. (Ill), but we also give results for additional 
strange quark masses on sets 1 and 2 and new results on sets 3, 6 and 8. We have 16,000 correlators per ensemble, except for 
sets 6 and 8 where we use approximately 10000. 



Set 








-a 3 (ipipi} 




aM Vg 


a fvs 




1 


0.013 


0.23637(15) 


0.11183(9) 


0.018607(29) 


0.0688 


0.53361(14) 


0.14199(6) 


0.045758(19) 












0.0641 


0.51491(14) 


0.13996(6) 


0.043616(19) 


2 


0.0064 


0.16615(7) 


0.10511(5) 


0.014524(18) 


0.0679 


0.52797(8) 


0.14026(3) 


0.045009(12) 












0.0636 


0.51078(8) 


0.13839(3) 


0.043038(12) 


3 


0.00235 


0.10172(5) 


0.09934(5) 


0.011762(11) 


0.0628 


0.50657(5) 


0.13720(3) 


0.042483(6) 


4 


0.01044 


0.19153(9) 


0.09075(5) 


0.011629(13) 


0.0522 


0.42351(9) 


0.11312(4) 


0.031756(10) 


5 


0.00507 


0.13413(5) 


0.08451(4) 


0.008511(9) 


0.0505 


0.41476(6) 


0.11119(2) 


0.030768(6) 


6 


0.00184 


0.08154(2) 


0.07988(2) 


0.006534(6) 


0.0507 


0.41481(2) 


0.11062(2) 


0.030768(4) 


7 


0.0074 


0.14070(9) 


0.06621(5) 


0.006153(8) 


0.0364 


0.30884(11) 


0.08238(4) 


0.019822(4) 


8 


0.0012 


0.05718(2) 


0.05781(3) 


0.002803(4) 


0.0360 


0.30483(4) 


0.08055(2) 


0.019504(2) 



configurations in each ensemble (somewhat fewer on sets 
6 and 8) and so the results have very small statistical er- 
rors. The valence light quark masses are equal to those in 
sea (except for a very small change on set 3) but we have 
shifted the valence strange quark masses slightly to be 
closer to the physical strange quark mass, following [2"5] . 
On sets 1 and 2 we give results for two different values 
of the strange quark mass, to help in constraining the 
valence mass dependence of the condensate. 

Errors on the condensate values are determined after 
binning over adjacent sets of at least five configurations, 
following analysis of the autocorrelation function. An 
example plot is shown, for coarse set 6, in Fig. [T] The 
autocorrelation function is defined as: 



Ca2 



(XjXi+AT) 



(x l )(x i+AT ) 



[.(■- 



(12) 



Here Xi represents a condensate value on configuration 
i and Xi+AT that on a configuration a further AT time 
units along in the ordered ensemble. AT = 1 thus cor- 
responds to adjacent configurations. Fig. [T] shows that 
nearby configurations in the ensembles are correlated and 
thus binning is necessary to obtain a reliable statistical 
error. A similar analysis applies to masses and decay 
constants as discussed in [55] . 

In the next section we describe the perturbative calcu- 
lation of the condensate which we will then subtract from 
the raw results of Table [TT] to enable the nonperturbative 
condensate to be determined. 



A. Perturbative calculation of (TrM -1 ) 

We computed the perturbative contribution (ifiip)pT to 
the chiral condensate for the HISQ action through first- 
order in a s : 

- a 3 (V>V>)PT,HiSQ = am x (13) 
[c (am ) + ci(am Q )a s + 0(a 2 s )] , 



where amo is the bare quark mass parameter that ap- 
pears in the HISQ action. The Feynman diagrams re- 
quired to this order are shown in Fig. [2] The perturba- 
tive quadratic ultraviolet divergence discussed in eq.([9]) 
shows up as finite values for the perturbative coefficients, 
as defined above, in the limit am — > 0. 




FIG. 2: Feynman diagrams for the calculation of the pertur- 
bative contribution to the quark condensate through 0(a s ). 



We computed the coefficients from numerical evalua- 
tion of the lattice loop integrals over a range of masses 
that includes the light and strange quark masses that we 
have used. A representative sample of our results is given 
in Table III and is illustrated in Fig. [3j 

An excellent fit to the perturbative coefficients in this 
range of small quark masses, aniQ < 0.1, can be obtained 
using the following parameterizations 



c (am ) = c 00 + {am Q ) 2 [c i log(am ) + c 02 ] , (14) 



and 



cx(am ) = cio + (am ) 2 [c n log 2 (am ) 
+ C12 log(amo) + ci 3 ] . 



(15) 



Higher order terms in am appear as discretisation errors 
in the comparison to MS to be done below and so can 
be ignored - they are negligible for the masses we are 
using in any case. The leading logarithm of amo a t each 
order originates entirely from the infrared region of the 



6 



loop momenta, and the respective coefficients coi and 
en can easily be computed analytically. The values of 
these coefficients must and do agree with the values in the 
MS scheme [115]. At one-loop we also have a constraint 
on the sub-leading (single) logarithm of amo since, as 
discussed in Appendix [Cj all log m terms must vanish 
in the difference between the vacuum expectation values 
of mi/iip in perturbation theory in the continuum and on 
the lattice. Allowing for the renormalisation between the 
A/5 mass and the HISQ bare mass: 



loga^ + 0.1143(3) 



m(/i) = m I 1 + 



(16) 

we find that c\i should have the value 0.2307(2). With 
the logarithmic terms fixed to their known values we can 
obtain the other coefficients in eqs. (TT4|) and ( f 5 ) from a 
fit to the values for co(amo) and ci(omo) as a function 
of am . We find: 



and 



coo 
coi 

C02 

ClO 
Cll 
Cl2 
Cl3 



0.38366(1), 

3/(2^ 2 ), 

-0.153(1), 

0.03657(7), 
-6/vr 3 , 
0.2307(2), 
0.308(15). 



(17) 



(18) 



These fits are illustrated in Fig. [3] and reproduce our 
results for the coefficients to within their numerical in- 
tegration errors, which are smaller than about 0.01% for 
Cq, and 1% for c\. 



The perturbative determination of the vacuum expec- 
tation value of tptp has also been done in the MS scheme, 
in |19j . The power divergence is missing in this case 
but, as discussed above, there are terms proportional to 
to 3 log m. ]T9 finds: 



[doil m + d 2 + a s {dull 
where l m = log(m(/z)//i) and 



(19) 



dl2lr, 



H3J 



doi 
do2 
du 
di2 
di3 



3 
2^2 



3 
4^2 



Cll 

5 



7T° 



5 
2^3 



(20) 



TABLE III: Zeroth- and first-order coefficients, Co and ci re- 
spectively, for the perturbative condensate, Eq. (14 1, for rep- 



resentative values of the bare quark mass parameter amo in 
lattice units. The uncertainties are from a numerical evalua- 
tion of the lattice perturbation theory loop integrals. 



amo 


CO 


Cl 


0.088 


0.37962(2) 


0.02561(12) 


0.079 


0.38029(1) 


0.02709(13) 


0.0728 


0.38075(1) 


0.02793(16) 


0.067 


0.38113(1) 


0.02902(18) 


0.062 


0.38146(1) 


0.02963(14) 




u.ooi I oyLj 




0.0505 


0.38212(1) 


0.03139(16) 


0.0448 


0.38240(1) 


0.03203(18) 


0.0386 


0.38271(1) 


0.03267(20) 


0.032 


0.38296(1) 


0.03349(24) 


0.028 


0.38313(1) 


0.03421(21) 


0.024 


0.38325(1) 


0.03424(21) 


0.020 


0.38336(1) 


0.03548(26) 


0.016 


0.38347(2) 


0.03545(34) 


0.01044 


0.38356(1) 


0.03614(49) 


0.00507 


0.38364(2) 


0.03631(31) 



As discussed in Appendix |B| we must subtract the dif- 
ference between the lattice QCD and MS perturbative 
calculations from our lattice QCD results to obtain the 
nonperturbative condensate in the MS scheme at the 
scale [l. We work with the combination mipi/j which 
would be RG-invariant in the absence of this perturbative 
contribution and it is convenient to derive the subtrac- 
tion needed in lattice units and as a function of the bare 



lattice quark mass. Using eq. (16) we obtain 



Apt = -a Um ipip) -p T ,mSQ - (m^V^PT.Sts, 

= c 00 (am ) 2 + a s c w (am ) 2 (21) 

+ (am ) 4 [c i^ - 0.077(1)] 

+ a s ( tt m„) 4 [c n ll + 0.1340(2)Z M + 0.406(15)] + . . . , 

where 1^ = log(/io). This difference of perturbative ex- 
pansions is now free of all logarithms of m and therefore 
well-defined and infrared safe. 



B. Determining the nonperturbative strange and 
light quark condensates 

1. A first look at the results 

The physical condensate in the MS scheme at the scale 
\x is then defined by: 



NP.MS 



( M ) = a- 4 (a 4 (m^) - A PT ) , (22) 



where (mifii/jja is the numerical result from lattice QCD 
and Apt is given through 0(a s ) as a function of the 
quark mass in Eq. (21). Apt will also contain unknown 
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y-axis in Fig. [4] is: 
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FIG. 3: Zeroth- and first-order coefficients, Co and ci respec- 
tively, for the perturbative condensate, Eq. (14 1, versus the 



bare quark mass parameter amo in lattice units. The un- 
certainties in Co resulting from numerical evaluations of the 
lattice loop integral are not visible in that plot. The fits given 
in the text are plotted as dashed lines. 



higher order pieces in a s that we can try to determine 
from a fit to the lattice QCD results. First we look at 
the effect of the calculated tree-level and one-loop con- 
tributions. 

Apt is a strong function of the quark mass in lat- 
tice units, dominated by the (ma) 2 terms that give rise 
to the quadrative divergence with inverse lattice spac- 
ing. This means that the relative size of the subtraction 
compared to the raw results varies strongly with quark 
mass and with lattice spacing, and this is reflected in the 
raw results before the subtraction is made. In Fig. [4] the 
open squares show the unsubtracted results (i.e. setting 
Apt to zero in Eq. (22)) as a function of the square of 



the inverse lattice spacing for quarks at the four differ- 
ent masses that we have results for in Table [TTJ strange 
quarks and light quarks of masses m s /5 (sets 1, 4 and 7), 
m s /10 (sets 2 and 5) and the physical value, m s /27 (sets 
3, 6 and 8). 

Instead of plotting the condensate results directly, the 



Ri 



for light quarks and 



R s 



Am s (t/jip s ) 



(23) 



(24) 



for strange quarks. The values of the raw unsubtracted 
R q are determined directly from Table |TT] using the am, 
aM, af and (ipip)o values given there. 

The ratio R q is a good quantity to plot (and later to 
use in our fits) for a number of reasons: 

• m(ipij>) is a physical renormalisation-group invari- 
ant quantity as m —¥ 0, up to discretisation errors, 
as is clear from the GMOR relation. The division 
by the square of the meson decay constant times 
its mass makes a dimcnsionless ratio which is con- 
venient but it is also one that (from the GMOR 
relation) we expect to be close to 1. 

• Using the ratio R q also reduces the effect of any 
slight mistuning of quark masses since the quark 
mass multiplied in the numerator cancels against 
the square of the meson mass in the denominator. 
The tuning of m s uses the ij s decay constant, as de- 
scribed in 25j. This means that, by definition, f Ve 
does not contain discretisation errors that would 
mask the identification of the pieces that diverge 
as a — > 0. 

• Finally the ratio has reduced finite-volume effects 
over that in either the numerator or denominator. 
This is expected from the fact that chiral loop ef- 
fects, which are sensitive to the volume, cancel in 
R q [27J This is illustrated in Fig. [5] where 
we show results [29] for pion mass, decay con- 
stant and (unsubtracted) light quark condensate as 
well as the ratio Ri of Eq. ( 23 ) for ensembles with 
am-Lsea = 0.00507 and am S:Sea = 0.0507 and three 
different spatial volumes. The spatial volumes cor- 
respond to a spatial length in lattice units of 24, 32 
and 40. The set with L/a — 32 is our set 4 (see 
Table [I]) . For each quantity we plot the ratio of the 
value at L/a to that at L/a — 40. It is clear from 
the plot that the finite volume dependence in each 
of m„, fir and (ipip)i is cancelled to a very high level 
of accuracy (0.1(1)% for set 4) in Ri. 

Figure [4] shows clearly the presence of a quadratic di- 
vergence with a~ 2 in the raw results. This is very 'clean' 
in our calculations because the form of the divergence 
is very constrained. Only a term of the form m q /a 2 is 
allowed in (rptp) for staggered quarks, i.e. no term of 
the form m 2 /a can appear. In the ratio R q this term 
takes the form Cm 2 /a 2 where C depends on the me- 
son mass and decay constant. The HISQ formalism has 
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FIG. 4: R q , defined as the ratio of quark mass times condensate in the MS scheme at 2 GeV to the square of the meson 
mass times decay constant, as a function of the square of the inverse lattice spacing. Left to right and top to bottom shows 
strange quarks and light quarks with masses m s /5, m s /W and the physical value. Squares use the unsubtracted condensate, 
pluses, the condensate after subtraction of the tree-level perturbative correction and crosses, the condensate after perturbative 
correction through one-loop. The value for ct s used to multiply the one-loop coefficient was a.y (2/a). Dashed lines illustrate 
very simple linear fits to the unsubtracted results as described in the text. 



very small discretisation errors, as is clear from the decay 
constant and meson mass results in |25j . and so there is 
little additional a-dependence to confuse the analysis of 
the divergent pieces. 

Because the power divergence is so dominant it is 
tempting to try to fit the unsubtracted results for R a 
to a very simple form: A + B/a 2 . This is in fact possible 
(it is important to include the error in the inverse lattice 
spacing when doing this since this is larger than the er- 
ror in Rq) and we obtain 1.02(3) + 0.725(3)/a 2 which is 
the dashed line in lefthand plot of Fig. |4j We also ob- 
tain 1.015(11) +0.229(5)/a 2 for Ri with m; = m s /5 and 
1.00(1) + 0.130(6)/a 2 for R t with m; = m s /10, shown 
in the next two plots in Fig. |4j These fits are too naive 
to be useful, as we shall see below, because they miss 
out many important terms. Consequently the value and 
error of the intercept, A, is unreliable for extracting a 
nonperturbative result for R q , especially in the s quark 
case. However, the fits do illustrate that the ratio of 



slopes is that expected for a term that behaves as m 2 /a 2 
(although the simple fit does not allow for the running 
of the lattice bare quark mass with scale). The ratio of 
slopes between that for s and for I with m; = m s /5 is 
3.2 which corresponds approximately to 5 (for the ratio 
of one power of the quark masses when the other power 
is cancelled by the square of the meson mass) divided 
by the ratio of the square of the decay constants from 
Table HH 

Figure H compares results for R q in which the tree-level 
piece of Apy has bee n su btracted from the raw values of 
m(V>V) following Eq. [22] We take n to be 2 GeV. These 
results are indicated by pluses. Now the slope in aT 2 is 
much smaller since the most of the divergence has been 
removed. This makes the results more sensitive to the 
form of the remaining pieces of the divergence and the 
simple linear fits that were made to the unsubtracted 
data are no longer possible. 

The 0{a. s ) perturbative contribution is very small for 
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FIG. 5: Finite volume effects in different quantities are il- 
lustrated by plotting the ratio of the quantity on lattices of 
spatial length, L/a, of 24 and 32 to that on lattices of spatial 
length 40. The lattices have the sea quark mass parameters 
of coarse set 5. The quantities shown are the pion mass (red 
pluses), pion decay constant (green crosses) and unsubtracted 
light quark condensate (blue bursts). Pink squares give the 
result for the quantity Ri defined in Eq. ( 23 1 [29] . Statistical 
errors (not shown) are approximately 0.1%. 



HISQ quarks and makes very little difference to the per- 
turbative subtraction. The crosses show the results tak- 
ing Apt to be the full calculated perturbative subtrac- 
tion through 0(a s 

*r 4 ' 



21 We have used 



given in Eq. 

at v ' *(2/a) |30j for the a s value multiplying the one- 
loop coefficient, but the coefficient itself is so small that 
variations in scale for a s make no difference. The crosses 
are barely distinguishable from the pluses giving the tree- 
level subtracted numbers. 

It is clear from Fig.[4]that there is still some divergence 
in a~ 2 left in R q after subtraction of the perturbative con- 
tribution through one-loop. This is not surprising since 
we know that Apt will have higher order terms in a s . 
The challenge now is to fit the one-loop subtracted R q 
allowing for these higher order terms and thereby obtain 
the physical, nonperturbative, results for the strange and 
light quark condensates. We will fit both R s and Ri si- 
multaneously and use the known mass dependence of the 
unknown higher order terms to constrain them. At the 
same time we will allow for higher-order non-divergent 
mass-dependent terms from perturbation theory as well 
as physical, non-perturbative, dependence on the quark 
mass. Possible dependence on positive powers of a, i.e. 
discretisation errors, must also be included. 



2. Determining a physical result from fitting 

We now describe the full fit to the results that we use 
to determine the final physical values for (ss) (= (ipip 3 )) 
and (11) (= (i/jipi)) a t the physical strange and light quark 
masses (where mi = (m u +md) /2). We take the following 



Rqfl{a, am q ) = i? N p phys + SR PT + SR a 2 + SR^ 



SRvoi- 
(25) 

Rq.o are the raw results obtained from Table [TT| i?NP,phys 
is the final physical result in the MS scheme at 2 GeV. 
The SR terms represent fitted or known dependence on a 
and am q . We use Bayesian techniques [31] to perform the 
fits so that we can add many higher order terms as part 
of each SR with constrained coefficients. This makes sure 
that the final error on i?NP,phys is not underestimated by 
ignoring the existence of higher order corrections. 

5Rpt contains the known tree- level and one-loop per- 
turbative results given in section |III A[ In addition we 
include unknown higher-order terms. For the aT 2 diver- 
gence these take the form: 



SR 



PT, div 



4<(«m ? ) 2 
ln {aUY{aM«) 



(26) 



with the analogous term for the strange quark case, with 
the same a n . a n is a coefficient whose prior we take to 
be 0.0 ± 4.0 and we allow for n = 2, 3 and 4. Note that 
a prior width of 4.0 is conservative given the size of the 
corresponding coefficient at tree-level and one-loop. For 
the non-divergent pieces we take: 



SR 



PT ,1 



4a!? (am q ) 4 



(27) 



again with the analogous term for the strange quark case, 
with the same coefficient c„. We take n = 2, 3 and 
4. c„ in principle contains a sum of powers of log(a/i) 
up to log n+1 (att). However, since log(a^) is small for 
fi = 2GeV and our range of lattice spacings, these pieces 
are negligible and do not affect the fit and we simply 
take a prior on c„ of 0.0(4.0). Again this is conservative 
given the results at tree-level and one-loop. For a s we use 
Q nf— (2/a) [3U] and discuss below the dependence of the 
results on changing 2/a to a different scale. a v f ~ 4 (2/a) 
takes values from 0.35 on the very coarse lattices to 0.26 
on the fine ensembles. 

SR a 2 allows for discretisation errors. We take the form 



Aa 



2i 



(28) 



Only even powers of a appear in discretisation errors for 
staggered quarks and we take their scale to be set by A w 
1 GeV. Since all tree level errors at 0(a 2 ) are removed 
in the HISQ formalism we take the prior for d\ to be 
0(a s ) i.e. 0.0(0.3). Higher order di are given the prior 
0.0(1.0). We include a 2 and a 4 terms, but have checked 
that higher order terms have very little effect. In addition 
we include a mass-dependent discretisation error in the 
form e(am) 2 , giving e a prior of 0.0(1.0). This allows for 
a number of effects, one of which could be mixing with a 
gluon condensate. This has negligible impact. 
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5R X includes the valence and sea quark mass depen- allow an additional linear dependence on the light quark 



dence that allows us to extrapolate to physical light quark 
masses and interpolate between the strange quark masses 
that we have to the physical strange quark mass. The 
chiral corrections to the GMOR relation were analysed 
in [27] (see also [28]). The leading corrections are par- 
ticularly simple because the chiral logarithms cancel to 
leave a correction proportional to Af 2 . We allow for both 
M 2 and M% terms in our light quark mass fits by defining 
a chiral expansion parameter 



Mi 



xi 



2(A X ) 2 ' 



(29) 



with A x = 1.0 GeV, and taking 



8IL. 



(30) 



We fit the m q — m s /5, m s /10 and rni tP hy S results with 
this form taking the prior on the gi coefficients to be 
0.0(2.0). This allows for a linear term of approximately 
the size expected in Higher order terms than xf 

have no effect. 

The chiral expansion of Eq. (30) combines with the 



-RNP,phys parameter in Eq. 25 to define the physical non- 
perturbative light condensate with its mass dependence. 
Since the GMOR relation is exact as m q — > we enforce 

this by taking the prior on R^p , to be 1.0000(5). Dif- 
ferences from 1 because of residual lattice finite volume 
effects up to 0.5% are allowed for as described in the 
paragraph above. 

The data for R s is fit simultaneously with that for the 
light quarks, because they share parameters for the per- 
turbative subtraction. However, we largely decouple the 
physical parameters because the strange quark is rela- 
tively far from the chiral limit and we would need a lot 
of parameters for a chiral expansion to connect light and 
strange quarks. Instead we allow a separate parameter 
for i4p. phys w ith the broad prior of 1.0(5). We take the 
same form for the valence mass dependence as in Eq. ( 30 ) 
where now 



(Ml - (0.6893(12)) 2 
2(A X )2 



(31) 



but this now simply allows for slight mistuning of the 
strange quark on some ensembles, and the fact that we 
have two values for the strange quark mass on sets 1 and 
2. 0.6893 is the physical value for the r] s mass deter- 
mined in [23]. The priors on gf' are the same as those 

on gf\ Finite volume effects are expected to be com- 
pletely negligible for R s because they are negligible for 
the components of R s . 

The strange quark results in Fig. [4] show some very 
small sensitivity to the sea light quark masses if we com- 
pare results on sets 1 and 2 and sets 3 and 4. We therefore 



mass in the sea in SR X of the form 



( Sm s 



V 10m 



s,phys 



(32) 



where 



(2m L 



+ m s 



(2m; ,phys + m s,phys )• (33) 



We take mi tP h ys = m s /27.4 ||2] and m StP h ys values de- 
termined from M Vs as in [23] . The prior for coefficient ki 
is taken as 0.0(1.0) which is conservative given the small 
effects observed in the results. We take SR x . sea to be 
common to both R s and R q for the light quarks. We 
note here also that the absence of chiral loop effects at 
this order means that staggered quark taste-changing ef- 
fects are also absent. They can be handled, if necessary, 
with a sea-quark mass dependent a 2 term [25] . Includ- 
ing such a term here makes no difference to the physical 
result. 

SR vo \ allows for remaining finite volume effects. These 
are small, as demonstrated in Fig. [3J They do produce a 
small systematic effect, however, because the lattice size 
in units of the pion mass, M n L, is somewhat smaller on 



the lattices with smallest m u / d . We take: 



fiRvol — 



— M L 



(34) 



where M is the pseudoscalar meson mass made of that 
quark (A/jr for Ri and M rjs for R s ) and L is the linear 
extent of the lattice from Table Q] Coefficient v is taken 
to have prior 0.0(0.2), consistent with Fig. [5j 

Fitting R s and R q , m q — m s /5, m s /l0 and rni tP h ys , 
results simultaneously to the form in Eq. ( |25[ ) readily 
produces good fits with x 2 /dof « 0.8 for 18 degrees of 
freedom. The final fitted result for R u /d and R s (eval- 
uated from R^ p phys and SR X taken at the appropriate 
physical masses) is robust to the addition of higher order 
terms in the various corrections. 

We take our final results from using 2/ a for the scale 
of a s . The results do not change significantly as this 
is varied (although the fitted coefficients a n do change). 
Our fits return a substantial value for the coefficient of 
the power divergent term at C(a 2 ), 02, of around 2.0, 
for d/a = 2/a. This is substantially larger than that 
seen at one-loop but not a particularly large value for 
a perturbative coefficient in general. It would simply 
imply that the small coefficient at one-loop for the HISQ 
action is not repeated at higher orders. We also find that 
the chiral correction to the GMOR for light quarks is 
substantial and negative = —1.7(6)). This will be 
discussed further below. 

The fit results are shown in Fig. [6] The data points 
(crosses) correspond to the lattice QCD results after sub- 
traction of the perturbative contribution through 0(a s ) 
(as in Fig. [4]). The filled bands show the fitted curves 
when the full fitted perturbative contribution (SRpt) is 
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TABLE IV: Error budget for the quantities R 3 ,phys, Ri.phys 
and their ratio defined in the text. Errors are given as per- 
centages of the final physical result. 





■R'S ,phys 


Rl,phys 


D 

- 1 L s ,pny s 


statistics 


6.1 


0.2 


5.1 


lattice spacing 


10.0 


0.3 


9.7 


finite volume 


1.5 


0.03 


1.5 


a s value 


1.7 


0.06 


1.7 


fitting power divergence 


7.5 


0.3 


7.2 


other perturbative subtraction 


1.3 


0.07 


1.3 


X&l extrap./interp. (m s ) 


3.0 


0.1 


2.9 


Xal extrap./interp. (mi) 


4.5 


0.2 


4.3 


a — > extrap. 


1.9 


0.05 


1.9 


sea mass effects 


0.5 


0.01 


0.5 


Total 


15 


0.5 


14.5 



subtracted and masses and decay constants are set to 
the physical values corresponding to the s quark and the 
light quark (for this we use M„ = M n o). These bands 
include the full error from the fit. 

Our final physical results for R q are the key results 
from this paper. 



Rl,phys 
Rs,phys 



4m ; (^^) ]g g(2GeV) 
4m,(^) ff g(2GeV) 

(fl M D 



(35) 



Wc find: 



-R>s,phys 
-R>l,phys 
Rs,phys 



R, 



l.phys 



0.574(86) 
0.985(5) 

0.583(84). 



(36) 



The complete error budget s for R s . p hys, Ri.phys and 
their ratio are given in Table IV The substantial 15% 
error that we have in R SjP hy S reflects the difficulty of ex- 
tracting a physical result from a power divergent quan- 
tity. For Ri the error is 17 times better largely because 
the slope of the divergent piece is 15 times smaller. Er- 
rors in R s , P hys are dominated by errors from the lattice 
spacing and from fitting the remaining power divergent 
subtraction terms. There are also substantial errors from 
statistics and from tuning to the light and strange phys- 
ical mass points. This is done by tuning the appropriate 
meson masses through the term SR x ^ va i in Eq. 30 This 



term depends on the lattice spacing through the defini- 
tion of xi (Eq. 29 1 and x s (Eq. 31), because the meson 



masses appear in GeV units in these terms. The uncer- 
tainties in these terms then becomes correlated with the 
fit to the power divergence, increasing the uncertainty. 
For Ri the power divergence is much less of an issue, but 
these same terms dominate the final error there as well. 
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FIG. 6: Results from fitting the ratio R for three different 
quark masses as described in the text. The crosses show the 
lattice QCD results after subtracting the perturbative values 
through 0(a s ). Black is for s quarks, blue for quarks with 
mass m 3 /5, red for quarks with mass m s /10 and green for 
quarks at the physical light quark mass. The dashed lines 
simply join the points of matching color for clarity. The filled 
bands show the physical curves for strange (black) and light 
(green) quarks, once the full subtraction of the fitted per- 
turbative contribution is made and masses are set to their 
physical values. The bands include the full error from the fit. 



The results for Ri and R s can be converted to values of 
the condensate using the lattice result for f Ve = 0.1819(5) 
GeV and M 7)b = 0.6893(12) GeV [25], and experimental 
values for / w (0.1304(2) GeV) and M w (0.13498 GeV). 
We obtain: 



m s (V^)f 6 (2GeV) 
roi{0V>)P*(2GeV) 



-2.26(34) x 10~ 3 GeV 4 
-7.63(4) x 10~ 5 GeV 4 . (37) 



The ratio of the two values above is slightly more accurate 
than a naive combination, giving 29.6(4.3). 

Using the precise determinations for light quark masses 
now available from lattice QCD we can finally obtain 
condensate values. We take mf s (2GeV) = 92.2(1.3) 
MeV [301 133] and mjmi = 27.41(23) [31 133- These 
give: 



(ss) MS (2GeV) = 
(ll) MB (2 GeV) = 



-0.0245(37) (3) GeV 3 
-(290(15) MeV) 3 

-0.0227(1)(4) GeV 3 
-(283(2) MeV) 3 , 



(38) 



where the second error for each condensate in GeV 3 
comes from the error in the quark masses. 

For the ratio of strange to light condensate we have: 



'(2 GeV) 



(70 MS (2GeV) 



1.08(16)(1), 



(39) 
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where the first error comes from R s / Ri and has the error 



budget given in Table IV and the second error comes 
from the strange to light quark mass ratio. 



3. Approach of R to the chiral limit 



The relationship of the light quark condensate to the 
chiral condensate is also important. R q is defined to have 
the value 1 from the GMOR relation in the chiral limit 
but the results of Eq. [36] indicate that it approaches this 
limit from below as the light quark mass is reduced. Sup- 
porting evidence for this is found by studying the quan- 
tity Rg derived from the combination of condensates used 
by the HOTQCD collaboration in their study of finite 
temperature QCD [35]. We define Rs by: 



Rs = 



Ami (Wl) 
f 2 M 2 ~ 



1 



(40) 



The quadratic divergence with lattice spacing cancels be- 
tween the two condensates because it is linear in the 
quark mass to all orders in perturbation theory. The 
non-divergent perturbative contributions proportional to 
the cube of the quark mass are completely negligible here, 
from the perturbative analysis in section III A| and so we 
do not need to include them in making Rs a physical 
quantity. Rs can then simply be calculated from the raw 
data in Table HU 
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FIG. 7: Rs as a function of m;/m s at three values of the lattice 
spacing. Points show the raw lattice results: the crosses are 
from very coarse lattices (sets 1 and 2 with two values of 
m s on each and set 3), open circles from coarse lattices (sets 
4, 5 and 6) and open triangles from fine sets 7 and 8. The 
shaded band gives the results of a simple fit incorporating 
discretisation and finite volume effects as described in the 
text. 



A plot of Rs against mi/m s is shown in Fig. [7] In the 
ttli — y limit on an infinite volume Rs — ► 1 as Ri does. 
Rs can be determined more precisely than Ri, however, 
because of the nonperturbative cancellation of the power 
divergence and it clearly approaches 1 from below. Rs 
differs from R[ by a term which is proportional to mi jm s 
and to the difference between (ss) and (11). Both the 
dependence of Ri on mj and the difference between Rs 
and Ri then contribute to the slope with mi seen in Fig. [7] 
We cannot separate them and therefore unambiguously 
identify the slope of Ri with m;. We can however use 
this for a consistency check. 

We fit Rs to the simple form: 



Rs = 1.000(1) +ci(l + c 2 a^+c 3 a 4 ) 



4 n m 



m s 



c 4 



mi 



c 5 e 



— m-n-h 



(41) 



This allows for linear and quadratic terms in mi with 
discretisation errors. We take priors on ci, C3 and C4 
to be 0.0(1.0) and prior on c 2 to be 0.0(0.3) consistent 
with a s a 2 behaviours. The final term allows for finite 
volume effects dependent on the combination m^L. As 
for our fit to Ri, we take the prior on C5 to be 0.0(0.2) 
for consistency with Fig. [5j 

The fit gives x 2 /dof of 0.75 for 10 degrees of freedom 
and a physical slope, c\, of -0.51(4). This value is consis- 
tent with the difference between Ri and 1 in Eq. |36| and 
indeed with the difference between R 8 and 1. This con- 
sistency between the results from Rs and Ri indicates 
that the difference between (ss) and (I/) (which would 
upset this consistency) cannot be large. This is indeed 
what we also find in Eq. [39] 



4- The chiral susceptibility 

A further quantity that is of particular interest in stud- 
ies of QCD at finite temperature is the chiral susceptibil- 
ity for a quark of flavor /: 



Xf 



d 
dm 



(-«%))• 



(42) 



We give results here for the chiral susceptibility for zero 
temperature QCD to fill out the physical picture of the 
condensate. From differentiation of the path integral for 
the condensate it is clear that the chiral susceptibility is 
given by the flavor-singlet scalar correlator. It is conve- 
nient to split this into two contributions which we call Xq 
and Xg ■ Xq comes from two scalar operators connected 
by quark lines, which is the flavor-nonsinglet scalar me- 
son correlator. x§ comes from two scalar operators con- 



1 In 1351 these are called Xconn and Xdisc- 
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TABLE V: Contributions to the chiral susceptibility, defined 
in Eq. (441 for very coarse set 1 and coarse set 4 for light 
(mi — m s /5) and strange quark masses. 



Set 


ma 


a 2 Xq 


a 2 Xs 


1 


0.013 


0.54296(36) 


0.045(14) 




0.0688 


0.45359(6) 


0.021(7) 


3 


0.01044 


0.50850(18) 


0.032(10) 




0.0522 


0.46231(3) 


0.014(4) 



nected only by gluons, in which the disconnected contri- 
bution is cancelled. 



X ~ Xq T Xg 



(43) 
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FIG. 9: Results for the gluon-connected contribution to the 
scalar correlator on set 1 for s quarks (blue crosses) and light 
quarks with m; = m s /5 (red bursts). The sum over time of 
this correlator is \g- 



The factors of number of tastes, N t , above are specific to 
naive/staggered quarks. 
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FIG. 8: Results for the quark-line connected scalar correlator 
on set 1: s quarks (blue crosses) and light quarks with mi = 
m s /5 (red bursts). The sum over time of this correlator is \q- 
Negative values of the correlator are not plotted on this log 
scale. 



Xq is readily calculated by generating quark propaga- 
tors with the same random wall source of noise as that 
used for the 7r and rj s mesons, but patterned with phases 
that are -1 on all odd sites on an even-odd partition- 
ing of the lattice. We then combine one of these prop- 
agators with the matching one used in the ir/r) s meson, 
again multiplying the odd sink sites with a phase of -1. 
Summing over spatial sites at each time slice gives the 
flavor-nonsinglet scalar correlator. Examples are shown 



in Fig. [8] Summing this correlator over time slices gives 
Xq- Results for \q f° r light and strange quarks on sets 1 
and 4 are given in Table |V| 

X g can be estimated from our existing results for the 
light and strange condensates. Because we have 16 time 
sources for our propagators we can determine the cor- 
relation between TrM -1 operator that are n time slices 
apart where n is a multiple of 3 for set 1 and a multiple 
of 4 for set 4. This gives a correlation function, for ex- 
ample that shown in Figure [9j Xg is then the sum over 
time slices of this correlation function. We can estimate 
Xg in several ways. Our central result comes from esti- 
mating an effective mass from the early time slices that 
dominate Xg an d where we have a strong signal. We can 
then reconstruct an estimated correlation function and 
sum over it. We can also simply sum over the correlator 
for the time slices that we have and multiply by 3 or 4 
as appropriate. From this range of methods we estimate 
the error in Xg as 30%. Values are given in Table |V| \g 
is much smaller than \ q . 



The quantities that are almost Renormlisation- Group 
invariant, that we can compare, are m 2 t Xf and 

mfi-itpijif). 



10 



m)xf 
Both of these 



This is done in Figure 
quantities contain the same power divergence (propor 



tional to m 2 f a 2 ) in the lattice spacing. In the suscepti- 



bility this divergence largely comes from Xg- The nondi- 
vergent perturbative contributions to the two quantities 
will be different but, as we have seen, they are very small 
and we ignore them here. 

We find that the difference between rnjxf and 

does not depend 



10 



—rrif^ipf), also plotted in Figure 
on the lattice spacing within errors. We simply average 
over the results at the two values of the lattice spacing 
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FIG. 10: Results for rnjxf (open circles) compared to 
mf{ip^)f) (open triangles) as a function of the square of the 
inverse lattice spacing, for s quarks (top) and light quarks 
with mi = m s /5 (lower plot). Results are for sets 1 and 4. 
Bursts show the difference of these two quantities. 



to obtain physical results for 



m f 1 - m/ 



d 



dirif 



-ipipf). 



(44) 



The values are: 2.18(7) x 10 3 GeV 4 for s quarks and 
3.93(9) x lCT 4 GeV 4 for light quarks with m ; = m s /5. 
Comparison of these num ber s with the physical results 
for mf(ipipf) given in Eq. (37 1 allows us to determine the 
value and sign of mjxf- For s quarks the comparison is 
straightforward and we find: 



0.08(35) x lO^GeV 



(45) 



consistent with zero. 

To obtain a value for mfxi a t m i — w s /5 we need a 
result for m^ipipi) at this mass. Our fit in Sec. IIIB2 
gives a result for Ri =s /§ of 0.915(26). At this value of mi 
we have M T = 315 MeV and can estimate /„• = 145 MeV 
from our results in Table [TTj Then mi—gj^ipipi—g/^) = 

4.78(20) x 10~ 4 GeV 4 . In the error we have included an 



interpolation error for each of the mass and decay con- 
stant of 1%. Then, subtracting the result for the differ- 



ence of Eq. 44 at m,[ = m s /5 given above, we have 



l l=s/5Xl=s/5 



0.85(22) x 10~ 4 GeV 4 , (46) 



which is a small positive slope. 

Our results are then consistent with a slope in (iptpf) 
with mf that is positive at m s /5 but may decrease or 
even become negative by m s . If the slope at m s /5 
remained constant for larger m/ it would give a total 
change in (ipipf) of 3 x 10~ 3 between m s /5 and m s which 
is not inconsistent with the change of 1.8(3.8) x 10 -3 that 
we find in Eq. |38| 



IV. DISCUSSION 

We have determined a physical value for the strange 
quark condensate from lattice QCD for the first time. 
This required both nonperturbative lattice QCD results 
and a perturbative determination of the power divergent 
contribution through 0(a s ). The calculation relies on 
the good chiral properties of staggered quarks to con- 
trol the form of the power divergence and the numerical 
speed and small discretisation of the Highly Improved 
Staggered quark formalism allow very precise results to 
be obtained at several values of the lattice with light u/d 
sea quark masses. We proceed by tracking deviations 
from the GMOR relation, making extensive use of our 
previous work determining the physical properties of the 
r] s meson, to obtain the strange quark condensate. Our 
best result comes from gluon field configurations that in- 
clude u, d, s and c quarks in the sea. The condensate is 
given in the MS scheme at a scale of 2 GeV. The evo- 
lution equation required to run m(ipip) to other scales, 
since it is not RG-invariant, is given in Appendix |Bj We 
obtain a very consistent result for the strange quark con- 
densate from independent calculations that include 2+1 
favors of sea quarks, as discussed in Appendix [D] 

Our value is —(290(15) MeV) 3 giving a ratio of strange 
to light condensates of 1.08(16). Earlier results come 
from sum rules of various kinds. These show signifi- 
cant variation and often have no estimate of the error 
associated with the value. Narison [35] gives a compila- 
tion with a final value for the ratio of strange to light 
condensates in the MS scheme at a scale of 2 GeV of 
0.75(12). A value of 0.74(3) is quoted from baryon mass 
splittings in |11] . Finite energy sum rules in the kaon sec- 
tor give a ratio 0.6(1) in |37| . More recently Maltman [T2"] 
uses sum rules for the ratio of decay constants JbJ Jb 
along with the 2007 lattice QCD average for this ratio 
of 1.21(4) [MIEH] to obtain (ss)/{ll} W g(2GeV) = 1.2(3). 
This updates an earlier result of 0.8(3) from Jamin [2"5] 
which used a quenched lattice QCD result for JeJ/b 
of 1.16(4). The current lattice QCD world average for 

f Bs /f B is 1.20(2) m- 
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FIG. 11: A comparison of results for the ratio of strange to 
light condensates in the MS scheme at 2 GeV. 



Fig. 11 compares our result for the strange to light con- 
densate ratio with the results from sum rules discussed 
above. Our central value lies between the sum rules re- 
sults, being in agreement with the larger value of |12j but 
only in marginal agreement with the lower values of [37] . 
A value below 0.6 is ruled out by our results at the So- 
level. Our value is more accurate than the result derived 
from fs 3 1 Jb and has the advantage over all the sum rules 
results that it is a direct determination from QCD and 
has a full error budget (Table IV). 



We obtain a very accurate value for the light quark 
condensate, giving —(283(2) MeV) 3 . We can distinguish 
the ratio Ri at the physical light quark mass from that of 
1 in the chiral limit using our results in Eq.( 36 ). Defining 
5„ from EH] 



Ri = 1 - S v 



(47) 



we obtain a value 5^ — 0.015(5). This is somewhat lower 
than the value of 0.047(17) estimated in [2S], although 
in agreement within 2a. It is not in agreement with the 
somewhat larger number of 0.06(1) obtained in [41]. Our 
result implies a value for the combination (2Lg — H^) of 
low energy constants from the chiral Lagrangian that is 
a factor of three lower than that used in [28] . The value 
is 3cr larger than zero, however. 

Note that we do not expect the value of the light quark 
condensate to agree with that of the chiral (zero quark 
mass) condensate, E. The relationship between them is: 



£(2 GeV) 



-(If) (2 GeV) f 2 (M 2 /m) rn=0 

Rl,phys f 2 M%/mi iphys 



(48) 



Here / is the decay constant in the zero quark mass limit 
and M 2 /m is the ratio of the square of the pion mass 
to the light quark mass in the same limit. f n /f can be 



determined from chiral extrapolation of lattice QCD re- 
sults. For example, a recent accurate calculation [3] gives 
f n /f = 1.0627(28) including 2+1 flavors of sea quarks. 
A preliminary analysis based on the results given here 
for 2+1+1 sea flavors in Table [jl| gives 1.056(1), in ac- 
ceptable agreement. From the figures in [3] we estimate 
(M 2 /m) m= o/(M 2 /mi tP h ys ) as 1.02. Combining these 
factors, along with Ri, into Eq. [48] makes clear that we ex- 
pect a 3% difference between the magnitudes of E 1 / 3 and 
((If)) 1 / 3 , dominated by the effect of f^/f (so that E is 
smaller). This is entirely consistent, assuming no differ- 
ence between 2+1 and 2+1+1 flavors of sea quarks, with 
the fact that we obtain (IZ)(2GeV) = -(283(2) MeV) 3 
and [3] obtain E(2GeV) = (272(2) MeV) 3 from a chiral 
analysis. 

Other methods for determining E are not as accu- 
rate, but in reasonable agreement. We quote here two 
recent examples. JLQCD/TWQCD give a result of 
(234(17) MeV) 3 [H] from the eigenvalue spectrum of 
overlap quarks with u, d and s quarks in the sea. The 
ETM collaboration give (299(38) MeV) 3 [43] from fits to 
the Landau gauge quark propagator with u and d quarks 
in the sea. Additional lattice results for E are collected 
in 03 . 

Our analysis has implications for other calculations. 
For example: 

• Finite temperature determinations of the chiral 
phase transition in QCD use an order parame- 
ter based on the light quark condensate. A non- 
perturbative subtraction is made with the aim of 
removing the power divergent pieces proportional 
to ma and with the assumption that the higher 
order ((ma) 3 ) terms are negligible. For example 
the HOTQCD collaboration uses an order param- 
eter [35] for visualising the transition (fits to find 
the transition temperature also include other re- 
sults) which is the ratio between non-zero and zero 
temperature of 



mi 
m s 



This quantity becomes 



mi 



l.NP 



,NP 



(49) 



(50) 



if we assume only the presence of a power diver- 
gence linear in ma. Our analysis shows that this 
is a good assumption. For example, the difference 
between subtracting only terms linear in ma at tree 
level and including terms cubic in ma is 0.2% for 
the strange condensate on the coarsest HOTQCD 
lattices (i.e. those with largest m s a values). An al- 
ternative might be to calculate (1 — md/dm)(ipip} 
as discussed in Sec. |IIIB4| The quark-line con- 
nected piece of this can be calculated directly by 



combining the expression for \ q in Eq. 44 and the 
expression for {iptp) in Eq. 10 The combination 
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becomes 

d_ 
dm 



"#)=2m £ TrlM^I 



(51) 



This can clearly be generalised to a sum over even 
source sites, implemented with a partial random 
wall. When combined with the quark-line discon- 
nected piece Xg this gives a physical quantity with- 
out power divergent pieces which is close to the 
value of the condensate itself. 

• The comparison of heavy-light current-current cor- 
relators to continuum QCD perturbation theory 
can be used to normalise heavy-light currents in 
lattice QCD. The light quark condensate appears 
in this comparison and the results given here will 
enable us to improve the analysis in [S]. This is 
underway. 

• Some recent papers [25] have speculated that the 
quark condensate may only be non-zero inside 
hadrons. A much smaller value outside hadrons 
would significantly ameliorate the fine tuning prob- 
lem associated with the cosmological constant. 
This suggestion appears to be in conflict with di- 
rect calculations of quark condensates as vacuum 
expectation values described here. 



V. CONCLUSIONS 

We give the first direct determination of the strange 
quark condensate from lattice QCD, having demon- 
strated how to extract a well-defined physical value from 
lattice results that contain a power divergence as the lat- 
tice spacing goes to zero. Our results include a calcula- 
tion through 0{a s ) in lattice QCD perturbation theory 
of the perturbative contribution to the condensate, part 
of which is the power divergence. The calculation relies 
on the good chiral properties of staggered quarks to con- 
trol the form of the power divergence and the numerical 
speed and small discretisation of the Highly Improved 
Staggered quark formalism to obtain precise results at 
multiple lattice spacings and light quark masses. Our 
results include values at physical light quark masses. 

We obtain a value for the strange quark condensate in 
the ES scheme at 2 GeV of -(290(15) McV) 3 . We give 
a full error budget for this result in Table |IV[ the main 
sources of error being those associated with fitting and 
subtracting the remaining power divergence. The result 
includes u, d, s and c quarks in the sea but we get good 
agreement with this value from independent calculations 
that include u, d and s sea quarks only. 

The value we obtain for the corresponding light 
quark condensate (where m; = (m u + md)/2) is 
-(283(2) MeV) 3 . Note that is significantly different from 
the value for the condensate in the chiral limit. The ra- 
tio of our light quark condensate to a recent lattice QCD 



value for the chiral condensate from jj] is 1.13(3), con- 
sistent with the behaviour of meson masses and decay 
constants approaching the chiral limit. 

We have shown that the ratio of four times the quark 
mass times condensate divided by the square of the me- 
son mass times decay constant approaches the GMOR 
value (of 1) from below as m; — > 0. At the physical light 
quark mass the value is 1.5% below 1, and at the strange 
mass it is 57% of 1. 

Our result for the ratio of the strange condensate to the 
light quark condensate is 1.08(16). This sits in the middle 
of the spread of results from QCD sum rules but provides 
significant additional information because it is a direct 
determination with a full error budget. The result will 
have impact on a number of other calculations both in 
the continuum and in lattice QCD. Some of the numerical 
techniques used here will be useful for determinations 
of, for example, the strangeness content of the pion or 
nucleon. 
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Appendix A: Condensates from correlators 



Eq. (Ill relates the zero- momentum pseudoscalar 



propagator to the scalar quark-condensate and is a well- 
known relationship |26j . The relationship is true (for the 
HISQ formalism) even on a single gluon configuration. 
Here, for completeness, we give a simple derivation, us- 
ing the equivalent naive quark formalism. 

The contribution to the propagator from a single gluon 
configuration is given by 



1 

7 + m 



1 



D-7- 



75 [ ^ ~ , | 75 

Ox \^ ' ' T "V xOJ 

(Al) 

where D is the gauge-covariant derivative, and the trace 
is over spin and color indices. The contribution to the 
scalar quark-condensate is given by 



-Tr 



1 



(A2) 



.^■7 + ^/00 
To extract the relationship, we multiply by the unit ma- 
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trix under the trace in the condensate: 



s = -J2 Tt 



■'•!./ 



(-D- 1 + m) 0x 
1 



1 

-D ■ 7 + m 



xy 



D ■ 7 + m 



.</<> 



(-^• 7 + m) 0;£ 



— (-D • 7) 2 + m / zo 



(A3) 



Only the m term in the numerator of the last expression 
survives the spinor trace since the other term results in 
traces of odd numbers of 7 matrices (which vanish) . Con- 
sequently 



S 



-mTr 



1 



X 



-(0. 7 )2 + m a /oii 
1 



1 



-D-j + mJ Qx \D-- 1 + mJ xQ 
1 



75 [ n r— I 75 x 

£>- 7 + mj 0x 



d -1 + ™J x q 



-mGps, 



(A4) 



which is Eq. (111. Since this relationship is true 



configuration-by-configuration, it must be true of the en- 
semble averages as well. 

Note that Eq. ( |A4| ) leads immediately to the GMOR 
relation (Eq. ([!])). To see this rewrite the pseudoscalar 
propagator in terms of its mesonic intermediate states. 
Only the pion contribution survives the m— >-0 limit, since 
the effective decay constants for excited states all vanish 
in that limit (by the Ward identity). The pion contribu- 
tion has an amplitude a — f 2 m%/(Am) multiplied by an 
exponential decay in time, whose integral gives 

The analysis of the propagator above only works for 
quark actions that have a 7^ piece and a scalar piece 
(and nothing else), and where those two pieces commute 
with each other. The commuting is essential if you want 



(-D ■ 7 + m)(D ■ 7 + m) = -(D ■ j) 2 + m 2 



(A5) 



with only terms having an even number of 7 M s on the 
right hand side. So this proof does not work for Wilson's 
lattice discretization of the quark action or similar formu- 
lations. On the other hand, it is true of staggered-quark 
formalisms such as HISQ. 



Eq.( 11 1 also follows directly from the (integrated) ax- 



ial Ward identity: 



£)<(m + m b ) J 5 ab (x)(m a + m b )J$(0)) = 

X 

-((m a + m b )(^a4>a + ^bipb)} (A6) 



with J^ b = ip a (x)j5tp b (x). This is exact on the lattice for 
lattice actions with sufficient chiral symmetry and again 
shows that is an identity, true configuration by configu- 
ration and for any m a and m b . See |46j for a derivation 
using twisted mass quarks. 

Here we have used the cases m a = m b and both equal 
to either mi or m s but we can derive from Eq.( A6) a 
relationship [33] between correlators for the mixed Gold- 
stone pseudoscalar made of light and strange quarks and 
the 'diagonal' cases: 

{ami + am s )^C K {t) = {ami)^C^{t) 
t t 

+ (am s )^C,,(t). (A7) 
t 

The left-hand side is then related to the sum of quark 
masses multiplied by the sum of quark condensates. This 
does not add new information so we do not make use of 
this relationship except as a test of our correlators. 

Appendix B: Condensates and the OPE 

Condensates typically arise in the non-leading terms of 
operator-product expansions (OPE). To illustrate, con- 
sider moments of two pseudoscalar densities composed of 
a heavy quark (mass M) and a light quark (mass m <C 
M) where the heavy quark fields are contracted with each 
other: 

{m + M) 2 { dxdtt n J 5 (x,i)J 5 (0) -> <D {n) (Bl) 



where 



O in) = / dxdtt n V>(x,i)75 



(m + M) 2 
D-j + M 



75^(0), (B2) 



and ijj is the light-quark field. Lattice simulations 
of (O|0< n >|O) can be used to determine the heavy quark's 
mass [47] . The (m + M) 2 factor makes (rl) independent 
of the ultraviolet regulator provided n>4; that is, lattice 
and continuum calculations should agree in the limit of 
zero lattice spacing. 

Operator is also short-distance, dominated by 

length scales of order 1/M, provided the heavy-quark 
is sufficiently heavy and the light quarks have momenta 
small compared with M . Consequently the OPE implies 
that can be expressed in terms of a set of local op- 
erators in an effective theory, with cutoff scale A < M, 
and coefficient functions that depend only upon physics 
between scales A and M: 



M n-4Q{n) 



1 (A) c[A/M,a s ,fn/M) 
(mi/iV) (A) 



iM 3 



d(A/M,a s ,fh/M) (B3) 
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where 1^ is the unit operator and we have replaced "0V 
by mipip, to simplify the coefficient function. Somewhat 
arbitrarily, we have chosen to express the right-hand side 
in terms of masses and couplings at scale /j, = M(fj,) = M: 

M = M(M), m = m{M), a s = a s {M). (B4) 



The effective theory on the right-hand side of Eq. (B3) 



could be, for example, lattice QCD with a lattice spacing 
a = 7r/A, or QCD with an MS regulator and fx = A. 

The coefficient functions c and d are perturbative 
when M is large, and analytic in a s and fh/M 2 . They 
can be computed using perturbative matching. For ex- 
ample, we can examine matrix elements of Eq. (B3) be- 
tween low-energy, on-shell light-quark states (q\ and \q'). 



The unit operator drops out and Eq. ( B3 ) can be rear- 
ranged to give 



d{A/M,a s ,fh/M) = 



rhM n - 1 {q\0 {n '>\q r 
(q\m^\q')( A ) 



PQCD 



(B5) 

where the right-hand side is computed order-by-order in 
perturbation theory. Since (q\Mipip\q') is independent 
of A, d is actually regulator independent: 



d = d(a s , m/M). 



(B6) 



Knowing d, one would then compute c using perturbative 
expansions of the vacuum expectation values: 



c(A/M,a s ,m/M) = 

/ <01QW |Q) (0|to^V|0) (A) 



V M 4 ~ 7 



fhM 3 



d(a s ,m/M) 



PQCD 

(B7) 



Eq. (B7| underscores the importance of avoiding 



normal-ordered operators in operator-product expan- 
sions. Each term on the right-hand side has infrared 
sensitive contributions that go like to 3 log (to). These 



cancel between the two terms in Eq. (B7|, order- by-order 



in perturbation theory; but this cancelation would have 
been ruined had we rep lace d ipip by the normal-ordered 
product : ipip : in Eq. (B3) (and c would no longer be 
perturbative) . 

It is also important to note that (O|to^-0|O)^ a ^ is not 
cutoff independent, because of operator mixing with the 
unit operator to 4 1, which implies that 



d(0|TO^|0) (A) 



dlog A 



7mix(a s (A),m(A)/A)m 4 (A). (B8) 



where 



7mix(Q: s ) 



3 
2^2 



Q(a s ). 



(B9) 



depends upon the regulator scheme beyond tree-level. 
This evolution is typically negligible for light quarks be- 



cause of the 



factor. 



Appendix C: MS condensates from the lattice 



The coefficient functions in operator-product expan- 



sions such as Eq. ( B3 1 are most conveniently computed 



using the MS regulator to define the operators on the 
right-hand side. On the other hand, the only technology 
available for determining the nonperturbative matrix ele- 
ments needed in such analyses is lattice simulation, using 
the lattice ultraviolet regulator. To combine these tech- 
niques we must be able to convert lattice determinations 
of (0|to?/>V'|0), for example, into the equivalent MS matrix 
elements. 

The relationship is again given by the operator product 
expansion: 



MS 



l (a)!! V/(^ Wa) 

+ ("#V)lqcd n / a ) 



(CI) 



where the coefficient functions / and h can only depend 
upon physics between fj, and the lattice cutoff n/a. In 
fact h — 1 since the matrix elements in 



(q\mi)tp\q r 



MS 



(g|m^v|g'>LQCD 



(C2) 



are fi and a independent, and therefore h must be a num- 
ber (and 1 is the correct number, from perturbation the- 
ory). The coefficient function / is computed order-by- 
order in perturbation theory using the expansions of the 
two condensates (computed with their respective regula- 
tors): 



/ = iL (< |m^|0)g- (0|m^|0)g ( 

n=0 

+ (aTO) 2 ^/W( aM)a =_( M ) 



PQCD 



(C3) 



n=0 



The cancellation of all log to terms between the two ma- 
trix elements in / is something of a miracle; it only works 
if the to in each case is precisely the to that mu ltiplies ipij) 
in the action for that case 3 . Apr in eq.( 21) is (am) 2 f 
calculated through 0(a s ). 



Additional terms appear in Eq. |C1| from mixing with 
higher dimension condensates, such as the gluon conden- 
sate. These are suppressed by positive powers of a. For 



2 We are ignoring nonperturbative short-distance physics (for ex- 
ample, small instantons) which can contribute to coefficient func- 3 If this isn't the case, then coefficient function h will not equal 
tions but is typically nonleading. one, but rather will be a series in 
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the gluon condensate the multiplier is (ma) 2 from chi- 
rality arguments. These terms then simply look like dis- 
cretisation errors in mipip and are handed as part of the 
general treatment of those errors. 



Appendix D: Lattice QCD calculation on nf 
gluon configurations 



2 + 1 



We show here further results for the strange quark 
condensate from two contrasting calculations, both using 
HISQ quarks, that include u, d and s quarks in the sea, 
but no c quarks. The first calculation uses sets of MILC 
configurations corresponding to 5 values of lattice spac- 
ing spanning a large range from 0.15 fm to 0.04 fm [48] 
and using the asqtad formalism for the sea quarks. The 
second uses HOTQCD configurations [35] and has more 
lattice spacing values (24 in total) but with only a lim- 
ited number (9) having accompanying meson masses and 
decay constants. The sea quarks are included using the 
HISQ formalism with u/d sea quark masses close to the 
physical value. The second calculation corresponds to 
the zero temperature results generated to accompany a fi- 
nite temperature analysis of the phase structure of QCD. 
This analysis needs many values of the lattice spacing for 
a fine-grained temperature scale, and the zero tempera- 
ture results are needed to fix the QCD parameters. The 
quark condensate is an important order parameter at fi- 
nite temperature but is also determined in [35] on the 
zero temperature ensembles. 

For the first calculation we use values of the strange 
quark condensate listed in Table [VT] These are obtained 
from studies of the rj s correlator on 9 different ensem- 
bles at 5 different values of the lattice spacing and mul- 
tiple sea quark mass values. The lattice spacing values 
we use here are defined from the rj s decay constant and 
are determined in 081 • From the values in Table VI we 



can construct the ratio R s defined in Eq. (24 1 and fit 



it as a function of lattice spaci ng in e xactly the same 

" The Q(a° ) and 



IIIB2 



III A 



way as that described in section 
0(a).) perturbative subtractions defined in section 
apply here also since no effects appear at this order from 
the differing number of sea quarks or the formalism used 
for them or the improvement coefficients in the gluon 
action (the MILC 2+1 asqtad configurations do not in- 
clude the rifa s a 2 improvement coefficients in the gluon 
action). These effects will cause differences in the per- 
turbation theory at 0(a 2 s ). The a 2 and higher order 
divergent pieces of the condensate are included in the fit 
with coefficients that, as before, are given a prior value of 
0(4). The appropriate a s value in this case is <Xy~ 3 (2/ a) 

rather than a^/~ 4 (2/a). Multiple valence s quark masses 
are given at each lattice spacing and we allow for lin- 
ear and quadratic dependence on the mistuning of the 
s quark mass, again as described in section |IIIB 2| We 
allow for mistuning of the sea quark masses through use 
of the parameter 5x sea [49] . 




0.03 



FIG. 12: Results from fitting the values for R s obtained 
on MILC configurations including 2+1 flavors of asqtad sea 
quarks, as described in the text. The crosses give the cal- 
culated values after perturbative subtraction through O(o s ). 
The hatched band corresponds to the fitted physical value 
after removing the remaining power divergence and discreti- 
sation and sea quark mass effects. Compare to Fig. [6] which 
includes 2+1+1 flavors of HISQ sea quarks. 



Figure [12] shows, as a hatched band, the physical re- 
sult from the fit, which has x 2 /dof = 0.4 for 20 degrees 
of freedom. For comparison the data points given are 
the values after perturbative subtraction through 0(a s ). 
The physical value obtained is 



Rs, P hy S = 0.555(84). 



(Dl) 



This is completely consistent with the result from 2+1+1 
flavors of HISQ sea quarks in section |IIIB 2[ and has 
a similar error. It is not such a complete calculation, 
lacking light quark mass results and not having such light 
sea quark masses, and is therefore not our preferred final 
result. It provides a strong check of our 2+1+1 result, 
however, being a completely independent set of numbers. 
The fits to the 2+1 results give very similar behaviour to 
that seen for the 2+1+1 case, for example choosing a 
coefficient of the a 2 /a 2 divergence of around 2. 

For the second calculation we use values of the strange 
condensate from the HOTQCD collaboration [55] , They 
generated ensembles with an improved gluon action and 
u/d and s quarks in the sea using the HISQ formalism. 
The QCD action differs slightly from that in section |TH] 
Apart from missing c quarks in the sea, the gauge field 
configurations here are improved through 0(a 2 ) at tree- 
level and without tadpole- improvement, i.e. the fairly 
substantial 0(a s a 2 ) improvement coefficients were not 
included. The lattice spacing was determined using the 
r\ heavy quark potential parameter, or tq on the coars- 
est lattices where r^/a < 2. The s quark mass was tuned 
by determining the mass of the r) s meson and the light 
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TABLE VI: Raw (unsubtracted) values for the strange quark condensate along with r] s masses and decay constants in lattice 
units calculated for valence masses given in column 4. The calculations use valence HISQ quarks on MILC configuration sets 
labelled in column 1 that include 2+1 flavors of asqtad quarks (see |48j for more details about the ensembles). The results 
are derived from the correlators calculated in |48| and [49] along with Eq.( 111, but we also give results for additional strange 
quark masses on sets 1 and 2. 8x sea is the mismatch between the sea 2m; + m s value and the physical result divided by the 
physical value of m s [49] . 



Set 


SGQ. 


a Vs (fm) 


Gilfh g Udl 


aM Vg 


af Vs 


-a 3 {tptp s } 


1 


0.47 


0.1583(13) 


0.061 


0.50490(36) 


0.1410(4) 


0.042399(38) 








0.066 


0.52524(36) 


0.1429(4) 


0.044637(38) 








0.080 


0.57828(34) 


0.1485(4) 


0.050795(37) 


2 


0.91 


0.1595(14) 


0.066 


0.52458(35) 


0.1434(3) 


0.044714(37) 


3 


0.64 


0.1247(10) 


0.0489 


0.41133(17) 


0.1124(2) 


0.030233(14) 








0.0537 


0.4310(4) 


0.1144(2) 


0.032423(15) 


4 


0.93 


0.1264(11) 


0.0492 


0.41436(23) 


0.1136(2) 


0.030585(21) 








0.0546 


0.43654(24) 


0.1160(3) 


0.033041(20) 








0.060 


0.45787(23) 


0.1182(4) 


0.035476(20) 


5 


1.5 


0.1263(11) 


0.0491 


0.41196(24) 


0.1135(2) 


0.030306(21) 








0.0525 


0.4259(6) 


0.1149(4) 


0.031817(23) 








0.0556 


0.4384(6) 


0.1161(4) 


0.033211(23) 


6 


0.59 


0.0878(7) 


0.0337 


0.29413(12) 


0.07954(9) 


0.018310(5) 








0.0358 


0.30332(12) 


0.08051(9) 


0.019273(5) 








0.0382 


0.31362(14) 


0.08171(15) 


0.020370(5) 


7 


1.1 


0.0884(7) 


0.0336 


0.29309(13) 


0.07959(11) 


0.018217(5) 








0.03635 


0.30513(20) 


0.08095(14) 


0.019467(7) 


8 


0.28 


0.0601(5) 


0.0228 


0.20621(19) 


0.0549(2) 


0.011311(5) 








0.024 


0.21196(13) 


0.0556(1) 


0.011851(3) 


9 


0.38 


0.0443(4) 


0.0161 


0.1525(2) 


0.0404(1) 


0.0075891(20) 



quark mass was taken as m s /20 (with some values avail- 
able for m s /5 but we have not used those). Results for 
the zero temperature strange condensate are available at 
24 values of the lattice spacing from 0.2 fm to 0.07fm 
(Table 14 of [35] gives values for two times the conden- 
sate). Note that these results were obtained by direct 
calculation of (TrM -1 ) using stochastic techniques. Cor- 
responding values of the lattice spacing are given in Table 
16 of [55!; some missing values can be inferred from the 
tables of temperature values at the corresponding value 
of /?. 



Figure 13 



shows the raw unsubtracted results for 
m a {ij}tl) s ) as a function of the square of the inverse lattice 
spacing, as well as the values after making the complete 
subtraction through 0(a s ) as given in Eq. (21). As for 
(Fig. [1]), the unsubtracted results 



i?, in section III B 1 



show clear evidence of a quadratic term in 



which 



is significantly reduced, but not completely absent, after 
the perturbative subtraction. 

For a subset of 9 lattice spacing values meson masses 
and decay constants are also given in Tables 18 and 19 
of [35] . In fact we use the 7 finest values only because the 
s quark mass is not as well-tuned for our purposes on the 
coarsest two lattices. Note that the decay constant values 
need to be multiplied by v2 to match the convention used 
here. For these we can construct the ratio R s given in [24] 
To obtain a physical result for R s we fit the subtracted 



> 
o 



A 
5- 




FIG. 13: Bursts give the raw data for m s {ipip s ) from [35] as 
a function of the square of the inverse lattice spacing. Open 
circles give results after subtraction of the 0(a a ) perturbative 
contribution. 



results as a f unction of lattice spacing in the same way 

3 (2/«) 



as m section 
rather than a 



IIIB2 



apart from the use of a?/ 



v ( 2 / fl )- 

The physical value for R s obtained from the fit is 



R s = 0.79(34) 



(D2) 



This is much less accurate than the result from sec- 



tion III but agrees both with that and the result from 
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the MILC 2+1 asqtad ensembles given earlier in this sec- 
tion. We have not extracted a light quark condensate 
from the HOTQCD results because finite volume sensi- 
tivity obscures the power divergence and leads to larger 
errors. 



We conclude from this that there is no sign of dis- 
agrecement between the strange quark condensate ex- 
tracted with u, d and s quarks in the sea and those in- 
cluding also c quarks in the sea. 
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